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Abstract 

By analyzing the potential energy surface, the shape phase diagram and the 
energy spectrum of the nucleus in U(5) symmetry in the IBM, we propose that the 
U(5) symmetry with parameters (A+B) < may be a model to describe the rotation 
driven vibrational to axially rotational phase transition along the yrast line. With 
such a model, we have described successfully the observed rotation driven shape 
phase transition along the yrast line of individual nucleus and proposed some other 
empirical evidences. 
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It has been well known that shape phase transition is one of the most significant topics 
in nuclear structure research. Many evidences of nuclear shape phase transition have 
been observed. For instance, in several isotopes, there exists shape phase transition from 
vibration to axial rotation or 7-unstable rotation with respect to the variation of neutron 
numberpp, and a triple point may appear [2 E|- Even in one mode of collective motion 
there may involve different characteristics, for example, along the yrast line there exists 
transition between rotations holding different relations between angular momentum and 
rotational frequency (referred to as band crossing and exhibiting a backbending) |3] . Very 
recently it was found that there involves rotation driven vibrational to axially rotational 
shape phase transition along the yrast line of individual nucleus jS]. On the theoretical 
side, the interacting boson model (IBM) has been shown to be successful in describing the 
shape phase transition along a series of isotopes[ri E]. And analytic solutions for the 
critical points of the phase transitions have been found [SI ElCni- The cranked shell model 
(CSM)PJ has been known to be able to describe the band crossing very well. However, a 
theoretical approach to describe the rotation driven shape phase transition from vibration 
to axial rotation along the yrast line in individual nucleus has not yet been established, 
even though several attempts have been made (see for example Ref.|l2j). By analyzing 
the potential energy surface and the energy spectrum of the U(5) symmetry in the IBM, 
we will show that the U(5) symmetry with a special choice of parameters can be a model 
to describe the rotation driven vibrational to axially rotational shape phase transition 
along the yrast line of individual nucleus. 

In the original version of IBM (IBMl), the collective motion of a nucleus is described 
by the coherent state of s- and rf-bosons, which hold angular momentum 0, 2, respectively. 
The corresponding symmetry group is U(6), and possesses three dynamical symmetry 
limits, namely, U(5), 0(6) and SU(3). With one- and two-body interactions among the 
bosons being taken into account, the Hamiltonian of the U(5) symmetry can be written 
asP, 

Hu{5) = Eo + 6dCiu{5) + AC2U{5) + BC20{b) + C'C'20(3) 5 (1) 
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where Ckc is the fc-rank Casimir operator of group G, and the parameters satisfy C ^ 
\B\ ^ \A\ ^ Ed- The IBM may be hnked to the collective model of Bohr and Mottelson by 
implementing the coherent state formalism|Tl ITBl IT^ I15j . In IBMl, the intrinsic coherent 
state of a nucleus with bosons is given by 



N 



|0) , (2) 



with 



ao = /?cos7, a±i = 0, a±2 = —^Psm'y. 
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It is evident that all components of (i-bosons can be described by the two parameters f3 and 
7, which have been shown proportional to the deformation parameters in the collective 
model pp. For instance, the relations (32 ~ 0.15/5 and 7 = 7 hold for rare earth nuclei, 
where P2 and 7 are the quadrupole deformation parameters in collective model. We then 
refer to the parameters (3 and 7 as deformation parameters in the follows. 

It is remarkable that the coherent state of Eq. (2) consists of all possible components 
of d bosons. As a consequence, it is not an eigenstate of angular momentum, but a 
superposition of states with all possible values of angular momentum. To investigate the 
shape phase structure and its transition among the states in ground-state band (or yrast 
band), we implement the technique of angular momentum projection |16tlT7| ITH] involving 
projection operator 



Pmk = / D^^^*{n)R{n)dn , (3) 



,L _2L + 1 

with M = K = 0. Then the energy functional [i.e. potential energy surface) of the state 
with angular momentum L in the ground-state band can be given as 

(A^;/3,7|Poo l^;/5,7) 
With the rotation operator written explicitly, Eq. (4) becomes 

p r o . £d/3'sm/3'd(^,i(3') {N;(3,^\He-^^''y [AT;/?, 7) 
^^'^ ' '^'^^ IodP'smP'dU(3'){N;(3,^\e-^^'^y\N;(3,^) ' 

where d^^, {(3') is the reduced rotation matrix. 



After some derivations, and considering that the deformation parameter (3 should be 
rather small for the vibrational to axially rotational (7 = 0) shape phase transition of 
current interest, we can approximate the potential energy surface of Eq. (5) up to (3^ as 

Egsb{N, L, (3) = Ao + \a{L - Lo)(3' + ^A^f3' + ^A^P' , (6) 



with 



Ao{N,L,ed,A,B,C) = ^[2ed + 8A + 6B + 4C + {A + B + 4C)L] L 

L(2N ~ L)(A + B) 
a = - 



2(3 + 2L) 
ed + 5A + AB 



-(A + B) ■ 

A4 and Aq are also functions of the interaction parameters Ed, A, B, C and the state 
indices A^, L. 

It follows from the spectrum generating principle that L G [0, 2A^], Sd > 0, while A, B 
and C may be either positive or negative. If {A-\- B) > 0, then a > and Lq < 0, so that 

— Lq) remains positive. Otherwise, if (A+B) < ( the dynamical symmetry condition 
guarantees {sd + 5^4 + 45) > ), then a < and Lq > 0, so that a{L — Lq) may change 
from positive to negative as the angular momentum L increases from under to over Lq. 
Because of the complexity of the A4 and Aq in terms of the parameters N, L, Ed, A, B and 
C, it is difficult to discuss the variation characteristics analytically. Numerical calculation 
indicates that, when the parameters in the Hamiltonian are taken as Sd > 0, {A + B) < 
and C > 0, the parameters A4, Aq in Eq. (6) can be negative, positive, respectively. It 
becomes then evident that the energy functional of Eq. (6) closely resembles the Landau 
free-energy ^211201 in the theory of thermodynamic phase transition, with L and f3 playing 
the roles of control parameter and order parameter, respectively. 

A more quantitative analysis is now in order for the potential energy surface of 
Eq. (6). If the parameters in the Hamiltonian (1) are chosen as (A + B) > 0, we have 
^ ^^'g^fji'^'^^ 1/3=0 = C({L — Lq) > 0. Then the potential energy is minimized only at /3 = 0. 
The nucleus with such a potential energy functional would always stay in a vibrational 
shape phase. If {A + B) < 0, then a < and Lq > 0, so ^-^^^WMl |/3=o > for L < Lq, 



and ^ ^3sb{N,L,i3) j^^^ < for L > Lq. The point /5 = evolves from a local minimum to 
a local maximum of the potential energy surface when L rises from below to above Lq. 
Furthermore, there exists a critical angular momentum Lc = Lq + j^r^, with which the 
energy surface acquires two and equal minima at /? = and P = ^SA^/AAq, respec- 
tively. There is another critical value L^ax for the angular momentum, beyond which a 
nucleus becomes unstable as the energy functional is no longer lower-bounded. It follows 
from Eq. (6) that Lmax = 2A^, the highest angular momentum possible, when Aq > 0, 
whereas L^ax = Lq + when A^ < 0. In any case, we have L^ax > Lq. For the states 
with L G (Lc, Lmax\i the energy functional is maximized at /5 = and minimized at some 
/9 7^ 0. As already mentioned, the energy functional has no minimum for the states with 
angular momentum L G {Lmaxi'^N]. The /5-dependence of the energy surface for some 
typical values of angular momentum L is illustrated in Fig. 1, where Lq and L^. are calcu- 
lated as 8 and 6 respectively. It becomes apparent that a nucleus has a stable vibrational 
shape phase with states around /3 = 0, if the Hamiltonian parameter {A + B) < and 
the angular momentum L < L^. In the regime of L G (LcL^ax], the only stable shape 
phase is rotational with states around the energy minimum at some /? 7^ 0. For L = L^, 
the energy surface has two degenerate but distinct minima, one of which is localized at 
(3 = 0, the other is localized at some (3^0. Under such a condition, the system may 
undergo a transition from one energy minimum to the other, and when that happens, the 
symmetry of the system is spontaneously broken. Such transition is of the first order, 
according to the standard theory of phase transition. The same theory would also predict 
the coexistence of vibrational and rotational shapes with L = as a precursor of a shape 
phase transition |2 H 1^. When the angular momentum L > Lmax, a nucleus may not be 
able to maintain a stable structure. Fig. 2 depicts the shape phase diagram of a nucleus 
in terms of the angular momentum L and the deformation parameter j3, where the same 
interaction parameters as in Fig. 1 are used, Lq and still are 8 and 6, respectively. 

The above discussion has shown explicitly that the potential energy functional derived 
from the U(5) symmetry with {A + B) < has the similar mathematical form of the 
Landau free-energy [121 120], which puts a U(5)-symmetric nuclear system in the standard 
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theoretical framework of first-order phase transition, so to correctly predict and well 
describe the vibrational to axially rotational nuclear shape phase transition in the ground- 
state band (or along the yrast line). 

On the phenomenological side, the energy spectrum of a nucleus in U(5) symmetry 
can be given as 

Eu(5) ^Eo + EdUd + A Udiud + 4) + S t(t + 3) + C L(L + 1) . (7) 

where n^, r, and L are the irreducible representations (IRREPs) of the group U(5), 0(5) 
and 0(3), respectively. 

Prom the spectrum generating process one can recognize that, if A > and B > 0, 

the states with r = Ud, L = 2nd form the ground state band and simultaneously the yrast 
band, and the energy of the state in the ground state band can be given as 

Egsbiud) = Eo + {A + B)nd^ + {sd + 4.A + 3B)nd + CL{L + 1) . 

The energy spectrum appears as the anharmonic vibrational one with increasing frequency 
huj = {sd + 4:A + 3B) + (A + B)nd (if the anharmonic effect induced by the rotation is 
taken into account, the vibrational frequency increases with Ud in the relation hcu — 
{£d + 4:A + 3B + 2C) + {A + B + AC)nd). However, if (A + S) < 0, the Eg,h{nd) {Eg,b{L)) is 
an upper-convex parabola against the Ud (L) and appears as the anharmonic vibrational 
one with decreasing frequency. Prom these anharmonic vibrational characteristics, one 
can recognize that, when (A + B) < 0, there exists a d-boson number n^^^ and an angular 
momentum 

L.^2„f'^-H(£l±^^-27V„. (8) 

where Nq = N with being the total boson number. As the angular momentum L > Lc, 
the yrast states are no longer the anharmonic vibrational ones mentioned above, but 
the quasi-rotational ones with Ud — Nq. Recalling the analysis in the coherent state 
formalism with angular momentum projection, one can conclude that the U(5) symmetry 
with parameters (A+B) < can describe the vibrational to axially rotational shape phase 
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transition along the yrast line. The Lc given in Eq. (8) is the critical angular momentum. 
For each yrast state with angular momentum L < L^, its energy can be given as 

Whereas for the one with angular momentum L > its energy should be expressed as 

E{L > L,) = E'o + CL{L + 1) , (10) 

and can be displayed as a part of an upper-concave parabola. For instance, for a system 
with = 15 and parameters Ed = 0.80 MeV, A = -0.025 MeV, B = -0.01 MeV and 
C = 0.004 MeV, with Eq. (8) we can fix the critical angular momentum Lc = 8h. The 
energy of the yrast states against the angular momentum L can be illustrated in the left 
panel of Fig. 3. It has been known that the energy of E2 transition 7-ray over spin (E- 
GOS) R = ■^t(-^"^-^~^) caji be taken as a quite good signature to manifest the vibrational 
to axially rotational phase transition along the yrast line[5^. As an auxiliary evidence, 
we show also the E-GOS of the yrast states with above parameters in the right panel 
of Fig. 3. The figure indicates apparently that the yrast states involve a vibrational to 
axially rotational phase transition and the angular momentum is definitely the critical 
point for the phase transition to take place. Such a transition is quite similar to that 
between the states with Up = and Up = A^|2^ in the vibron model with random 
interactions. 

Reviewing above analysis, one knows that, in the anharmonic vibrational model, there 
involves competition between vibration and rotation. In the case of that the interaction 
parameters are taken as {A + B) < 0, the vibrational frequency decreases and the rota- 
tional effect increases if the angular momentum (or d boson number) increases. As the 
angular momentum L (or Ud) reaches the critical value, the vibration disappears, so that 
only the rotational effect governs the property of state. Simultaneously, a sudden increase 
happens in the d boson number. Then the structure of the wavefunction changes from the 
vibrational one to the rotational one, so that the vibrational to the axially rotational shape 
phase transition takes place. In the microscopic point of view, cranking random phase 



approximation calculation has shown that the gradual decrease of vibration can induce 
a backbending, i.e., a shape phase transition[25 . On phenomenological level, the U(5) 
symmetry with parameter A < has been used to describe the collective backbending of 
high spin states successfully^]. 

As an application of our presently proposed model, we analyze the typical exam- 
ple ^°^Ru involving the rotation driven vibrational to axially rotational shape phase 
transition 5 . By fitting the experimental data of the yrast band of ^°^Ru with Eqs.(9) 
and (10), we obtain the energy spectrum and the E-GOS plots as shown in Fig. 4. To 
show the sudden change of the wavefunction, we also list the configuration [n^, r, L) of 
the sates in the yrast band in Fig. 4 and the fitted parameters in Table 1. The figure and 
the Table show evidently that the shape phase transition happens at angular momentum 
L = 12 and our model describes well such a shape phase transition. 

Table 1: Fitted parameters of the bands ^°^Ru, ^^^Cd, ^^^Cd, "^Te, ^^^Sm and ^^^Hg 



band 


Bd (MeV) 


A (MeV) 


B (MeV) 


C (MeV) 


N 


Lc 




0.5281 


-0.009376 


-0.008819 


0.01475 


20 


12 


ii^cd 


0.6581 


-0.007640 


-0.01327 


0.01440 


25 


10 


ii^Cd 


0.5695 


-0.0003789 


-0.01586 


0.01393 


29 


10 




0.7747 


-0.01678 


-0.01423 


0.01269 


17 


12 


i42Sm 


1.101 


-0.06263 


-0.01165 


0.01713 


7 


10 


188 Hg 


0.5957 


-0.03057 


-0.002976 


0.009193 


9 


12 



We have also analyzed the available experimental energy spectra of the yrast bands 
of even-even nuclei with 30 < Z < 100 and simulated the data with least-square fitting 
in our present model. We find that, besides the yrast bands of nuclei ^°^Ru, ^^^Cd and 
^^^Cd identified by Regan and collaborators in Ref.|5j, the yrast bands of nuclei ^^^Te, 
^^^Sm and ^^^Hg involve the vibrational to axially rotational shape phase transition. The 
theoretical results of the energy spectra and the E-GOS plots of these bands and the 
comparison with experimental data are illustrated in Fig. 5. The fitted parameters are 
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listed in Table 1. From Fig. 5 and Table 1, one can infer that our model can describe the 
rotation driven shape phase transition successfully. 

Looking over the fitted boson number listed in Table 1, one may know that, for ^^^Sm 
and ^^^Hg, it agrees well with the simple ansatz: boson number is half of the valence 
nucleons (or holes). However, for the nuclei in A ~ 110 mass region, it differs from 
half of the valence nucleons obviously. Recalling the recent suggestion that the valence 
neutrons may be in the /iii/2 orbitaljSl 123; we infer that the microscopic configuration 
of the nuclei in A ~ 110 mass region is so complicated that the boson number of them 
can not be simply taken as half of the valence nucleons. It means that the effective 
boson number may be important to the structure of these nuclei. In fact, quite early 
IBM calculations j2Hl 12^1 ISHl HI] had shown that, to describe the spectroscopic property 
of these nuclei well, effective boson number should be implemented. Then we believe 
that the fitted boson number is reasonable, even though it is not consistent with half 
of the valence nucleons for the A ~ 110 nuclei. Meanwhile it provides a clue that the 
structure of the nuclei involving rotation driven vibrational to axially rotational shape 
phase transition is very complicated and needs sophisticated investigation. 

In summary, by analyzing the potential energy surface of the nucleus in U(5) symmetry 
in the coherent state formalism with angular momentum projection in the IBM, we give a 
phase diagram of the vibration and the rotation in terms of the angular momentum and 
the deformation parameter in this letter. We have then proposed that the U(5) symmetry 
with parameters {A + B) < may be a model to describe the rotation driven vibrational 
to axially rotational shape phase transition along the yrast line. With such a model, we 
have described successfully the vibrational to axially rotational phase transition along the 
yrast line in A ~ 110 mass region identified by Regan and collaborators [S]. By analyzing 
the available experimental spectra of even-even nuclei with 30 < Z < 100, we show that, 
besides the ones in A ~ 110 mass region, the yrast band of nuclei ^^^Te, ^^^Sm and ^^^Hg 
may also be the empirical evidences involving the rotation driven vibrational to axially 
rotational shape phase transition. 
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Figure 1: The energy surface (Landau free energy) of a nucleus against the "deformation 
parameter" j3 at some typical angular momentum L (with the parameters in Eq. (1) 
being taken as sa = 1.010812 MeV, A = -0.153669 MeV, B = 0.0822402 MeV, C = 
0.0185084 MeV. ). 
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Figure 2: The vibration and rotation phase diagram of a nucleus in terms of the angular 
momentum L and the "deformation parameter" /3 (with the same parameters for Fig. 1). 
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Figure 3: An example of the energy against spin (left panel) and the E-GOS (right 
panel) along the yrast line (filled squares and triangles, respectively) in the approach of 
U(5) symmetry with (A + B) < (with parameters £d = 0.80 MeV, A = -0.025 MeV, 
B = —0.01 MeV, C = 0.004 MeV. The solid and dashed lines are implemented to guide 
the eye.) 
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Expt. 


Theor. io2pj^ 
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Figure 4: The theoretically obtained energy spectrum and E-GOS plot of the yrast band 
of ^''^Ru(curve) and the comparison with experimental data (filled squares, taken from 
Ref.[Sl) 
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Figure 5: The theoretically obtained energy spectrum and E-GOS plot of the yrast bands 
of ^^^Cd, ^^^Cd, ^^^Te, ^^^Sm and ^^^Hg (curve) and the comparison with experimental 
data (filled squares, taken from Refs. |^ IHH IM| ) 
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